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We consider a graphene sheet folded in an arbitrary geometry, compact or with nanotube-like 
open boundaries. In the continuous limit, the Hamiltonian takes the form of the Dirac operator, 
which provides a good description of the low energy spectrum of the lattice system. We derive 
an index theorem that relates the zero energy modes of the graphene sheet with the topology 
of the lattice. The result coincides with analytical and numerical studies for the known cases of 
fuUerene molecules and carbon nanotubes and it extend to more complicated molecules. Potential 
applications to topological quantum computation are discussed. 
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Introduction:- Much attention has been focused lately 
on various geometric configurations of graphene, where 
an interplay takes place between geometry and electronic 
properties such as its conductivity P, [3, H, 0, H, S 01 ■ P'^s- 
vious methods for obtaining the zero modes of the sys- 
tem (electronic eigenstates with zero energy) are based 
on lengthy analytical or numerical procedures. As a pos- 
sible alternative the much celebrated index theorem Q 
offers an analytic tool that relates the zero modes of ellip- 
tic operators with the geometry of the manifold on which 
these operators are defined. This theorem has a dramatic 
impact on theoretical and applied sciences 0. It allows 
to gain information about the spectrum of widely used 
elliptic operators by simple geometric considerations that 
could be otherwise hard or even impossible to determine. 

It is the purpose of the present letter to establish a 
version of the index theorem that relates the number of 
zero modes of graphene wrapped on arbitrary compact 
surfaces to the topology of the surface. Nanotube-like 
open boundaries that are of relevance to physical config- 
urations of graphene are also presented. When consid- 
ering the low energy limit of graphene a linearization of 
the energy is possible due to the presence of individual 
Fermi points in the spectrum. This results in a Dirac 
equation defined on the manifold of the lattice [3| , which 
describes the low energy behavior of the system well. An 
additional coupling to an effective gauge field with long 
range effect is generated by the deformations of the lat- 
tice needed to introduce curvature. Since the Dirac op- 
erator is an elliptic operator, it is possible to employ the 
index theorem 0, [lo| to obtain information about the 
low energy behavior of graphene and in particular about 
its conductivity properties. Indeed, as we shall see in 
the following an exact relation can be found that con- 
nects the number of zero modes with the genus, <?, and 
the number N of possible open ends of the surface. Our 
results are in agreement with the known cases of icosahe- 
dral fuUerene molecules [HI and graphite nanotubes 
where the spectrum has been determined analytically or 



numerically. A relation between the zero modes of more 
complicated molecules is provided. 

There is a variety of applications that spring from 
this work. Information about the spectrum of complex 
molecules constructed out of nanotubes can be provided 
that may be impossible to obtain with other analytical 
approaches. Moreover, the presence of G fermionic zero 
modes dictates the existence of a 2^ ground state de- 
generacy of the initial Hamiltonian. Hence, one could 
employ reverse engineering and construct a fermionic lat- 
tice Hamiltonian with a particular degeneracy structure. 
This is of much interest in the area of topological quan- 
tum computation 1^, 14 1, where information can be en- 
coded in the degenerate states protected by topological 
considerations. Moreover, lattice Hamiltonians such as 
the one that models graphene can be engineered by Fermi 
atom gases superposed with optical lattices . The de- 
tection procedure of the zero modes of these models is 
well established flG*! offering an alternative method for 
probing the conductivity properties of various graphene 
configurations. Similar approaches for the ground state 
degeneracy of fractional quantum Hall systems in the pla- 
nar case or on hi gh-genu s Riemannian surfaces have been 
taken in [B, [H, [liM HH . 

The model:- Let us first consider graphene, a flat iso- 
lated sheet of graphite. It can be shown [4] that the 
tight-binding approximation reduces the system to that 
of coupled fermions on a honeycomb lattice (see Fig.[T]). 
The relevant Hamiltonian is given by 
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<i,j> denotes nearest neighbors on the 
ai are the fermionic creation and annihi- 



where i > 0, 
lattice and a| 

lation operators at site i with the non-zero anticommuta- 
tion relation {ai, a]} = <5y ■ The corresponding dispersion 
relation can be easily evaluated as 



E{p) = ±iWl -^4cos^ 
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where the interatomic distance is normaUzed to one. As 
can be deduced from ([2]) at half-filhng, graphene pos- 
sesses two independent Fermi points instead of Fermi 
hnes. This rather unique property makes it possible to 
hnearize its energy by expanding it near the conical sin- 
gularities of the Fermi points. It is not hard to show that 
the resulting Hamiltonian is given by the Dirac operator 
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where the Dirac matrices, 7", are given by the Pauli 
matrices, 7" = cr", and ± corresponds to the two inde- 
pendent and oppositely positioned Fermi points. Hence, 
the low energy limit of graphene is described by a free 
fermion theory. The corresponding spinors are given by 
(I K±A) , I K±i?))^ where A and B denote the two differ- 
ent sublattices that comprise the honeycomb lattice (see 
Fig. [T]) and Kj- denote two independent Fermi points 
chosen such that K_ = — K+. One can see from ^ 
that any fermionic mode rotated by is mapped onto 
another mode with the same energy and opposite mo- 
mentum. This fact does not necessarily hold for the zero 
modes of non-flat geometries as we shall see in the fol- 
lowing. 




FIG. 1; The honeycomb lattice comprises of two triangular 
lattices, A, denoted by black circles and, B, denoted by blank 
circles. A single pentagonal deformation can be introduced by 
cutting a 7r/3 sector and gluing the opposite sites together. 



Curvature deformations:- In order to evaluate the ef- 
fect of curvature on the fermionic modes it is instruc- 
tive to consider first the transformation properties of the 
spinors. A 27r/3 rotation of the lattice centered on a 
hexagon is given by R2-k/3 = exp(i^(T^), which is an 
SU(2) element acting on the spinors, while a tt rota- 
tion Rtt = exp(i^(T^)zT^ mixes both spinor elements and 
Fermi points, where is the Pauli operator acting on 
the K± components [H, H, The form of i?^ is due 
to a reorientation of the reference frame to its original 
direction; it distributes an i and a —i to the spinor com- 
ponents while the rotation of the Fermi momenta takes 
an additional minus sign indicated by ir*'. 

To obtain surfaces with arbitrary topology, curvature 
is introduced to an initially flat honeycomb lattice by in- 
serting deformations. In doing so, we shall demand that 



each lattice site has exactly three neighbors and that 
the lattice is inextensional (free to bend, impossible to 
stretch). The minimal alteration of the honeycomb lat- 
tice that can introduce curvature without destroying the 
cardinality of the sites is the insertion of a pentagon or 
a heptagon; this corresponds to locally inserting positive 
or negative curvature, respectively. Other geometries are 
also possible, leading to similar results. 

To introduce a single pentagon in a honeycomb lat- 
tice, one can cut a 7r/3 sector and glue the opposite sides 
together, as illustrated in Fig. [1] This causes no other 
defects in the lattice structure. We shall demand that 
the spinors are smooth along the cut remedied by intro- 
ducingcompensating fields which negate the discontinu- 
ity [1,11] . Indeed, the cut introduced in Fig. [T] causes an 
exchange between A and B sublattices. This can be in- 
corporated into the Hamiltonian by introducing the non- 
abelian gauge field. A, with circulation 



A" 2 



that mixes the -f and — spinor components. This flux 
can be attributed to a fictitious magnetic monopole in- 
side the surface with a charge contribution of 1/8 for 
each pentagon [22] ■ In addition, moving a frame around 
the deformation gives a non-trivial transformation that 
is equivalent to a spin connection Q. The flux of this 
field around the pentagon is given by 
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and measures the angular deficit of 7r/3 around the cone. 
These fields exactly compensate the spinor transforma- 



tion -Rtt/s = ^TrR2Tr/3 — exp(— «^cr^)ir'' produced when 
the lattice is rotated by tt/S. 

The modified Dirac equation, which incorporates the 
curvature and the effective gauge field, couples the K± 
spinor components together due to the non-abelian char- 
acter of A. They can be decoupled by a single rotation 
that gives 
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where fc = 1 , 2 denotes the components in the rotated 
basis with ^ A'^dx^ = ±7r/2. is the zweibein of the 
curved surface with metric g^^i/ that defines the local flat 
reference frame, ijafs = e^e^g^i/. This equation faithfully 
describes the low energy behavior of graphene, such as its 
zero modes, when it is deformed to an arbitrary surface. 

Index theorem:- We have seen how our system reduces 
to the Dirac equation of a spinor field on the surface of a 
lattice interacting with a gauge field. Our aim now is to 
construct an index theorem that gives a relation between 
the zero modes and the particular topology of the surface 
on which the graphene sheet is wrapped. 
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Let us briefly review the index theorem. Consider a 
compact, oriented, smooth, Riemannian manifold M and 
the elliptic operator D over M. Here D is the Dirac oper- 
ator given in ^ . One can show that for compact M the 
Dirac operator is self-adjoint with a discrete spectrum 
of eigenvalues [23|. For even dimensional M, such as a 



surface, the spinor space breaks up into two irreducible 
pieces, denoted by and V~. The Dirac operator in- 
terchanges between the two spaces in the following way 
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Denoting by i'± the dimension of the zero eigenspace of 
V^, one can define the index of D by 

indcx(Z?) = — 

One can introduce the operator 75 = i Y[a 7" with eigen- 
values ±1 that breaks the Hilbert space into the two sub- 
spaces ® V" with the property tr(75) = dim(F+) — 
dim(T^~) = — While there is a spectral symmetry 
between the non-zero modes, the nuU-subspace does not 
need to be symmetric due to topological defects, a fact 
that makes the index(£') a non-trivial quantity. Further- 
more, the index theorem @, S, M, ^ states that 



mdex(D) — — 



F, 



(6) 



where F = d /\ A is the field strength of possible gauge 
interactions and the integration is taken over the whole 
surface. Note that there is no curvature contribution to 
the index theorem for two dimensional manifolds. 

Our aim is to evaluate the contribution from the gauge 
field, F, in It is solely determined from the geomet- 
ric characteristics of the surface through the presence of 
pentagons and heptagons on the original lattice. We shall 
only consider a surface that is either compact or when 
open boundaries are present then a smooth differentiable 
surface can be produced when the surface is "glued" with 
its mirror symmetric one. One can calculate the number 
of deformations in a lattice necessary to generate such 
a surface by employing the Euler characteristic. Indeed, 
for V, E and F being respectively the number of vertices, 
edges and faces of a lattice on a surface with genus g, 
and TV open ends the Euler characteristic is given by 

X = V - E + F = 2(1 - 5) - iV. 

We can easily verify that a single cut in the surface can 
reduce the genus of the surface by one and increase the 
number of open ends by two , i.e. {g, N) (,g — 1, iV+2), 
thus preserving the Euler characteristic, x- From the 
total number of pentagons, hexagons and heptagons in 
the lattice, n^, rtg and ny, respectively, we see that E 
= (5n5 + 6ne + 7n7)/2, V= (5^5 -t- Gng -I- 7ri7)/3 and 
F= + riQ + rij, giving finally 



This reflects the fact that pentagons and heptagons have 
equal, but opposite curvature and gauge fiux contribu- 
tions, while non-trivial topologies necessarily introduce 
an imbalance in their numbers. 

Eqn. ([7| reproduces the known case of a sphere (x = 2, 
n5 = 12 and riy = for the Ceo fuUerene), a torus (x = Oi 
"-5 = "-7 = for the nanotubes) or the genus-2 (x = —2, 
71-5 = Oj '^7 = 12) where equal numbers of pentagons and 
heptagons can be inserted without changing the topology 
of the surface. 

Now we are in position to evaluate the index(_D). The 
contribution of the gauge field term in ^ can be cal- 
culated straightaway from the Euler characteristic. It is 
obtained by adding up the contributions from the sur- 
plus of pentagons or heptagons. Thus, the total flux of 
the effective gauge field can be evaluated by employing 
Stokes's theorem, giving 
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clef. 
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where ± corresponds to the fc = 1, 2 gauge fields. Hence, 
from one obtains 



3x/2, for fc = 1 
-3x/2, for fc = 2 



(8) 



As both of the cases contribute zero modes to the system 
the least number of zero modes is given by 3|x| — |6— 6(7— 
3N\, which coincides with their exact number if = 
or = 0. 

This result reproduces the number of zero modes for 
the known molecules. The fuUerene, for which genus g ~ 
and = has six zero modes which correspond to the 
two triplets of Cgo and of similar larger molecules [i , 24 1 . 
For the case of nanotubes we have g = and N = 2, 
which due to formula ^ gives v+ — v- = 0. This is 
in agreement with previous theoretical and experimental 
results 
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ng - ny = 6x = 12(1 ~ 5) - 6iV. 



(7) 



An alternative derivation of index(D) for the case of 
nanotubes is given when considering the operator 75, 
which in our case is given by 75 = ia^a^ = —a^. As we 
have seen, for the case of flat graphene sheets, transfor- 
mations with respect to a" give modes with exactly the 
same energy even for the null subspace. When nanotubes 
are considered then the dispersion relation ([2]) is accom- 
panied by the boundary condition SN^Px + VSNyPy = 
27rm, where Nx, Ny determine the relative lattice po- 
sition of the sites that are identified, and m is an in- 
teger [2^. One can easily see that the 75 symmetry 
between the zero modes is still preserved giving, as ex- 
pected, zero difference between them. This symmetry 
also holds when additional boundary conditions are em- 
ployed to generate, e.g. the torus. In the case of the 
fuUerene, the topological defects cause the breakdown of 
this symmetry. 
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Conclusions:- The presence or absence of zero modes 
in physical systems is of much wider interest. For exam- 



ple, spin lattice models are proposed [26j, |27[ that exhibit 
ground state degeneracy that is unaffected by small per- 
turbations and, thus, capable of supporting error free 
quantum information encoding. In our case G fermionic 
zero modes contribute 2*^ ground state degeneracy, where 
G depends on the topology of the surface of the lattice. 
As we have seen, minimal local deformations that do not 
change the topology can be introduced by having equal 
numbers of pentagons and heptagons added in the lat- 
tice. These deformations do not alter the number of zero 
modes given by the index theorem. Hence, this method- 
ology presents a promising way to construct topological 
models with protected ground state degeneracy. 

Finally, we would like to present an alternative phys- 
ical model that can realize Hamiltonian ([1]). It consists 
of a single species ultra-cold Fermi gas superposed with 
optical lattices in a honeycomb lattice obtained by three 
planar standing wave lasers [28j . Recent experiments [l5l | 
demonstrate that it is possible to control this system to a 
high degree of accuracy obtaining very low temperatures 
of the order of Q.VTp for arbitrary filling factors, where 
Tp is the Fermi temperature. At half filling one can re- 
alize the dynamics of llj near the Fermi points, thus, 
simulating the conducting properties of graphene in pla- 
nar configuration. At this point one could systematically 
study, for example, the effect of disorder. The latter can 
be introduced cither by deforming the geometry of the 
lattice, e.g. by introducing pentagonal configurations at 
the boarders of the lattice, or by considering impurities 
and lattice defects |29|. 

Detection of zero modes in similar systems has alr eady 
been achieved in the laboratory in the following way 
Consider the setup where the Fermi gas is trapped with 
the optical lattice and an additional wide harmonic po- 
tential. When the system is brought out of its equilib- 
rium, i.e. translated from the minimum, Xmin, of the 
harmonic trap, then two distinctive cases can arise de- 
pending on the presence or absence of zero modes. The 
conducting regime is characterized by cloud oscillations 
around Xmim while in the insulating case the center of 
mass of the cloud remains at the displaced position ex- 
hibiting small Bloch oscillations. This approach can of- 
fer an alternative experimental verification of the relation 
between topological defects and the conductivity proper- 
ties of graphene. 
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